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ABSTRACT 

We consider the approximation of optimal discrete-time linear quadratic Gaussian (LQG) 
compensators for distributed parameter control systems with boundary input and unbounded 
measurement. Our approach applies to a wide range of problems that can be formulated in a state 
space on which both the discrete-time input and output operators are continuous. Approximating 
compensators are obtained via application of the LQG theory and associated approximation results 
for infinite dimensional discrete-time control systems with bounded input and output. Numerical 
results for spline and modal based approximation schemes used to compute optimal compensators 
for a one dimensional heat equation with either Neumann or Dirichlet boundary control and 
pointwise measurement of temperature are presented and discussed. 
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1. Introduction 


In this paper we develop an approximation theory for the computation of optimal discrete-time 
linear quadratic Gaussian (LQG) compensators (combined feedback control law and state estimator) 
for distributed parameter systems with boundary input or control and unbounded output or 
measurement. In a continuous time setting, boundary input typically results in an unbounded input 
operator. That is, the system's input operator maps the control input into a space larger than the 
state space in which the open-loop system is usually formulated. In the discrete-time case, on the 
other hand, for a wide class of distributed systems, the resulting input operator is bounded on the 
usual underlying state space. By unbounded output or measurement is meant that the system 
output operator has domain in a space smaller than the usual open-loop state space. 

For continuous time systems, Pritchard and Salamon (1987) have established an abstract 
semigroup theoretic framework for treating the linear quadratic regulator problem (control only) for 
infinite dimensional systems with unbounded input and output operators. Their approach is based 
upon a weak or distributional formulation of the Riccati equations which characterize the optimal 
feedback control laws in an appropriate dual space . Curtain (1984) provides a procedure for the 
design of finite dimensional compensators for parabolic systems with unbounded control and 
observation. In (Curtain and Salamon, 1986) a finite dimensional compensator design procedure 
for a wider class of infinite dimensional systems with unbounded input (but bounded output) 
including hereditary systems with control delays and partial differential systems with boundary 
control is developed. Lasiecka and Triggiani have looked at linear regulator problems for parabolic 
(1983a, 1987a) and hyperbolic (1983b, 1986) systems with boundary control and obtained, 
among other things, global and local regularity results for the optimal controls and state 
trajectories. In (Lasiecka and Traggiani, 1987b) Galerkin approximations and an associated 
convergence theory for closed-loop solutions to regulator problems for parabolic systems with 
Dirichlet boundary input are studied. A more complete survey of the boundary control literature 
including references to some of the poineering work in this area can be found in (Pritchard and 
Salamon, 1987). 
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In our treatment here, we consider the discrete-time problem (i.e. piecewise constant input and 
sampled output). Our interest in the discrete-time or digital formulation is motivated by 1) the fact 
that it represents a more accurate or realistic description of how the linear-quadratic theory for 
distributed systems would actually be applied in practice, and by 2) how the boundedness of the 
discrete-time input operator in the usual underlying state space facilitates the development of an 
approximation theory which can simultaneously handle both unbounded input and unbounded 
output. Our approach is based upon an application of the theory we developed earlier in (Gibson 
and Rosen, 1985 and 1986) for the approximation of optimal discrete-time LQG compensators for 
infinite dimensional systems with bounded input and output. Our results are applicable to a rather 
wide class of boundary control systems in which a restriction of the state transition operator and the 
discrete-time input operator are bounded on a space on which the output operator is bounded as 
well. 

An outline of the remainder of the paper is as follows. In Section 2 we describe an abstract 
framework for the study of boundary control systems and their discrete-time formulation. In 
Section 3 we review the LQG theory for infinite dimensional discrete-time systems and associated 
abstract approximation results. In the fourth section, we illustrate the application of our technique 
on an example involving the spline and modal subspace based approximation of optimal 
compensators for a one dimensional heat equation with either Neumann or Dirichlet boundary 
control and pointwise measurement of temperature. Section 5 contains some concluding remarks. 


2. The Boundary Control System and its Discrete-Time Formulation 

We employ a semigroup theoretic formulation for a class of abstract boundary control systems 
which has appeared elsewhere in the literature. See, for example, (Curtain and Salamon, 1986). 
Let W,V and H be Hilbert spaces with W and V densely and continuously embedded in H. We 
consider boundary control systems of the form 
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(2.1) 

w(t) = Aw(t), 

t>0 

(2.2) 

O 

£ 

II 

/—N 

o 


(2.3) 

Tw(t) = v(t). 

t>0 

(2.4) 

y(t) = Cw(t), 

t>0 


where A e 33(W,H), the boundary input operator T is an element in £ (W,R m ) and the output 
operator C is an element in 33(V,RP). 

We require the assumptions that 1) T is surjective and its null space, TU(T) = {tp e W: 

Tcp = 0}, is dense in H, 2) the operator <U, defined to be the restriction of the operator A to Tl (T). 
is a closed operator on H and has non-empty resolvent set and 3) for each T > 0, all w 0 £ W, and 
v e C^O.T; R m ) with Twq = v(0), there exists a unique w e C([0,T]; W) fl (^([O.T]; H) which 
depends continuously on Wq and v and which satisfies (2.1) - (2.3) for each t e [0,T]. It then 
follows (see Hille and Phillips, 1957) that the operator 0, : Dom (Cl) c H — > H given by Cl<p = 
A(p for tp eDom(Cl) = 71 (T) is the infinitesimal generator of a 13 0 semigroup, (<T(t): t > 0}, 
of bounded linear operators on H. 

Define the space Z as the dual of Dom(Cl )c H. Then H is densely and continuously 
embedded in Z, { <T (t) : t £ 0} can be uniquely extended to a 13 q semigroup of bounded linear 
operators on Z, and the operator Cl can be uniquely extended to the operator in 33(H,Z) given by 
( Cl <p)(\j f) =«p, Ci*\j/> H for cpeH, \|reDom(Cl*). 

For v e R m , since T is onto, there exists a w e W such that Tw = v. With this w define 
Bv = Aw - Clw. It can be shown (see Curtain and Salamon, 1986) that Hi is a well defined 
element in 33 (R m , Z) and that for each w 0 e H and v e L 2 (0,T; R m ) there exists a unique w e 
C([0,T]; H) D H 1 (0,T; Z) which depends continuously on w 0 and v and which satisfies 

w(t) = Clw(t) + Bv(t), t > 0 
w(0) = w 0 

in Z. The function w is given by 

t 

(2.5) w(t) =ST(t)w 0 + J ^(t - s) Bv(s)ds, t>0 

o 
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and is referred to as a weak solution to the boundary control system (2.1) - (2.3). 

The discrete- time formulation of (2.1) - (2.4) is found by considering piecewise-constant 
controls of the form 

(2.6) v(t) = u k , t e [kx, (k + 1 )t), k = 0,1,2,... 

where x denotes the length of the sampling interval. Letw k = w(kx), k= 0,1,2,... where w( ) 
is the unique weak solution to (2.1) - (2.3) given by (2.5) corresponding to w 0 e H and input v 
given by (2.6). (We note that with piecewise constant input of the form (2.6), the solution w is in 
fact a strong solution on each subinterval [kx, (k+l)x].) Recalling that T is a surjection, let 
e £(R m , W) denote any right inverse of T and for each k = 0,1,2,... define 
^ e C([kx, (k + l)x]; H) by z^t) = w(t) - T+u k , t e [kx, (k + l)x]. Then 
z k (t) = w(t) = Ciw(t) + Bu k 
= Ci ^(t) + (a + BDr+u k 
= Clz^t) + A Pujj , 1 6 (kx, (k + l)x], 

Z k (kx)= W k - Fu k . 


Therefore 


w k + i =z k(( k+1 ) x ) + r+ u k 

= 0" (x)(w k - p-u k ) + (s) A r f u k ds + r+u k 

= ^T(x)w k + (I - 2T(x)) r f u k + f 0" (s)AT H 'u k ds, 


or 


w k+l = Tw k + Bu k’ k = 0,1,2,... 

w 0 eH 
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where T e 25(H) and B e £(R m , H) are given by T = U'(t) and B = (I - O' ft)) r+ + 
rT 

J <T (s)Ar + ds respectively. 

We note that the discrete-time input operator B is well defined and does not depend upon a 

particular choice for P\ Indeed if Bj and B 2 are the input operators which correspond to the 

choices r + and T + then for ueR m we have 
1 2 

(Bj - b 2 )u = a - <r (x»( r* - r+)u + f tr (s)A (r+ - r+>uds. 

But (r+ - T+ )u e TUT) = Dom(Ct) and therefore 

X X 

J tT(s)A(r+ - n )uds = | 3'(s)Ct(p--r )uds 
0 0 

X 

= Ji Sr(s)(r r _i 2 )uds= (5r(x> _i)(i t _r 2 )u - 
0 

In addition, if r* - is chosen so that c 'll, (A), B takes on the particularly simple form B = 

(I - <T (t))r f . It is also worth noting that a simple calculation reveals that 

B = Jq $T(s) Bds 

in agreement with the standard technique for obtaining the discrete or sampled time formulation of 
a continuous time system in either a finite dimensional or bounded input setting. 

It is our intention here to apply the approximation theory we developed earlier in (Gibson and 
Rosen, 1986) for the design of optimal discrete-time LQG compensators for infinite dimensional 
systems with bounded input and output operators. We therefore require the additional assumptions 
that 4) T = <T(x) e £(V) and 5) RCr 1 ") c V. Although not all boundary control systems we 
might formulate would satisfy these conditions, as will become evident in Section 4 below, a wide 
class of interesting and important systems do. In this case, the control system (2.1) - (2.4) takes 
the form 
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(2.7) 

( 2 . 8 ) 
(2.9) 


w k+1 = Tw k + Bu k , k = 0,1,2,... 
w 0 e y 

y k = Cw k , k =0,1,2,... . 


3. LOG Theory for Infinite Dimensional Discrete-Time Systems and Finite 
Dimensional Approximation 

The discrete-time linear quadratic control or regulator problem for the boundary control system 

(2.1) - (2.3) takes the form 

Find u* = (u k ) k _ 0 e i 2 (0> °°> R m ) which minimizes the quadratic performance index 

J(U)= 2 < Q w k- w k > v + u k Ru k 
k=0 

where Q e £(V) is self-adjoint and nonnegative, R e £(R m ) is a symmetric positive definite 
mxm matrix and the state w = { w k) k _ 0 ev °l ves according to the recurrence (2.7), (2.8). 

The optimal control is given in closed-loop, linear state feedback form by u* = -Fw*, 
k = 0,1,2,... where F = (R + B IIB) ^^ T and n is the minimal nonnegative, self-adjoint 
solution (if it exists) to the operator algebraic Riccati equation 

(3.1) n = T*(n - ITB(R + B*IlB)' 1 B*n)T + Q . 

A control u is said to be admissible for the initial data w 0 if the resulting state trajectory w k = 
w k (w 0 ,u), k = 0,1,2,... is such that J(u) < oo . If for each w 0 e V there exists an admissible 
control u, then the Riccati equation (3.1) admits a self-adjoint nonnegative solution IT. If, in 
addition, an admissible control drives the state w k to zero, asymptotically as k — > «>, then this 
solution is unique (see (Gibson and Rosen, 1985)). The optimal state trajectory w* = { w k } * 
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evolves according to w* = S k w 0 , k = 0,1,2,... where the closed loop state transition operator 
S e 2 j(V) is given by S = T - BF. If Q is coercive, then S has spectral radius less than one and is 
uniformly exponentially stable. From the finite dimensionality of the control space we obtain 

(3.2) u* = - <f,w*> v , k = 0,1,2,... 

k k v 

m 

where f = (f 1 ,f 2 ,...,f m ) T e x V is referred to as the optimal functional feedback control gains. 

j=i 

When only a finite dimensional measurement y = {y k } Q of the infinite dimensional state w 
is available (i.e. equation (2.9)) a state estimator or observer is required. For a given input 
sequence u and corresponding output sequence y the optimal LQG estimator is given by 

(3.3) w k+1 = Tw k + Bu k +F{y k -Cw k }, k = 0,1,2,... 

(3.4) w 0 e V 


where the optimal estimator or observer gains F e £ (R p ,V) are given by F = TIIC (R + 
CIIC )'* with II e £(V) the minimal self-adjoint, nonnegative solution (if one exists) to the 
operator algebraic Riccati equation 

(3.5) fi = T( fi - n c*(R + c n cVc n yr* + q . 


A A 

The operator Q e £(V) is assumed to be self-adjoint, nonnegative and the pxp matrix R is 

A 

assumed to be symmetric, positive definite. In a stochastic setting, the operator Q and the matrix 

A 

R are assumed to be respectively the covariance operator and matrix for uncorrelated, zero-mean, 
stationary, Gaussian white noise processes which corrupt the state and measurement. In this case, 

A 

if Q is trace class, (3.3), (3.4) is the infinite dimensional analog of the discrete-time Kalman-Bucy 

A A 

filter. In a strictly deterministic setting, Q and R are assumed to be determined via engineering 
design criteria such as stability margins, robustness of the closed-loop system, etc. 
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Sufficient conditions for existence and uniqueness of self-adjoint, nonnegative solutions II to 
the algebraic Riccati equation (3.5) are of course analogous to the ones given earlier for the control 
problem (i.e. with regard to the usual duality which exists between the optimal LQG control and 
estimator problems). 

A 

Since F e (R p ,V), it has a representation of the form 
Fy = f T y , y e R p 

* T P 

where f = (f r f 2 ,..., f J 1 e x V is referred to as the optimal functional observer gains. 

j=i 

When the optimal control law is used together with the optimal observer, that is 
(3.6) u* = -F w* , k = 0,1,2,... 

where w = { )^° 0 is given by (3.3), (3.4) with input u = u and corresponding output 
y = y*, we obtain the optimal LQG compensator. The resulting closed-loop system is given by 

<Uf k = /6 k Tir 0 , k = 0,1,2,... 

where *Uf k = (w k , w* ) T with {w k } k = 0 the state trajectory which results from the input (3.6) 
and A> e £(VxV) is given by 

't -bf 

A = 

FC T-BF-FC 
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If we define e k to be the difference between w k and w* ; e k = w k - w k > we find e k = S^q > 

A /V 

k = 0,1,2,... where S = T - FC. It can be shown that if S and S are uniformly exponentially 

A 

stable, so too is A and its spectrum , g(A), is given by a(A) = o(S) U o(S ). 

For each N = 1,2,... let V N denote a finite dimensional subspace of V and let P N be a 
bounded linear mapping from V onto V N (for example, the orthogonal projection with respect to 
either the V or H inner product). Let T N , Q N , Q N e 3S(V N ), B N e£(R m ,V N ) and C N e 
£(V N ,R p ) and set 

f n = (R+B*n N B N )'B*n N T N 


and 


*vl 


f - t n n N c N (R + c N n N c N > 


where II N and II N are the minimal, self-adjoint, nonnegative solutions (once again if they exist) to 
the finite dimensional operator algebraic Riccati equations 

(3.7) n N = T-( n N - n N B N (R + B^n N B N )-> b^tt n )t n + q n 


and 


(3.8) 



A <{{ 

respectively. The approximating optimal compensator is given by u* = -F N w_ k = 0,1,2,.. 

v% JN 

where w * = { w * } 00 is determined according to the approximating observer 

N N,k k = 0 
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w" 


Nk+1 


= T N w* + B n u* + F N {y* - C N w* }, k = 0,1,2.. 

N Nk ^ MV 1 ' V M I, IN MV 


N k 


H k 


N k 


W 


N,0 


= P 


N w 0 


e V 


N 


The measurements 


* 



are given by 



Cw N>k , k 


0,1,2,... where 


w N >k ^= Tw N,k- Bu * k - k = 0,1,2,... 
w N,o = w o- 


The resulting closed-loop system is given by W N k = ^ N ^ n, 0’ ^ = where 

1^N,k = ( w N,k- w * k ) T and 6 ^ (VxV N ) is given by 


(3.9) 


= 


-BF, 


N 


FC T -B F • f C 
r N^ N D N r N 


Let " T n ' B n F N and “ t n ‘ Fn^n and assume that the spaces V N are V-approximating 
in the sense that P N ->I strongly on V as N — > <» . Assume further that 
T n P n -> T, t;p n -» T*, Q n P n -4 Q and Q N P N Q strongly on V and that 
B n — » B and C N P N — > C in norm as N — » «*>. If the pairs (T N , B N ) and (T^ , C* ) are 

A 

uniformly exponentially stabilizable and the pairs (T N , Qn) and (T* , Q N ) are detectable (see 
Kwakemaak and Sivan, 1972) then there exist unique, self-adjoint, nonnegative solutions n N and 

A A 

II N to the algebraic Riccati equations (3.1) and (3.5). If II N and II N are bounded from above 

A A 

uniformly in N, then II N P N and n N P N converge weakly to II and II respectively as N — » <*> . 

A 

If, in addition, S N and S N are uniformly exponentially stable, uniformly with respect to N, 

A 

then II N P N and II N P N converge strongly. Weak convergence of n N P N to II yields strong 
convergence of F N P N to F and S N P N to S. If II N P N converges strongly then F N Pn —> F in 

A A A A A 

norm. Weak convergence of II N P N to II yields weak convergence of F N to F and S N P N to 
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AAA A A A A 

S. When n N P N — » IT strongly, then F N — > F in norm and S N P N — » S strongly in V as 
N — » Finally, if we let <P N denote the mapping of VxV onto VxV N given by tP N (Wj, w 2 ) = 

(Wj,P N w 2 ) then II N P N -» II weakly or strongly is sufficient to conclude that 

A A 

weakly or strongly depending only upon whether n N P N — > n weakly or strongly as N — > <» . 
Under appropriate additional hypotheses on the spectral properties of the open-loop system, the 
nature of the approximation spaces v N and the mappings P N’ it is possible to obtain a result 
regarding the norm convergence of /6 N <P N to A. Norm convergence of the closed-loop state 
transition operators is sufficient to conclude that uniform exponential stability of A implies uniform 
exponential stability of /6 N for all N sufficiently large (see Gibson and Rosen 1986). 

In practice, the finite dimensional approximating subspaces V N are often constructed using any 
one of a number of common finite element bases, e.g. polynomial and hermite spline functions, 
mode shapes, orthogonal polynomials, etc. For some of the discrete-time boundary control 
systems of particular interest to us here, the approximations to T and B, t n and B n , are typically 
obtained by approximating the continuous time semigroup, { <T(t) : t > 0}, by a semigroup of 
bounded linear operators on V N , { <T N (t): t > 0}. (In actual fact it is the infinitesimal generator Ct 
of the semigroup { <T (t): t > 0} which is approximated by a bounded linear operator Ct N on v N with 
{ *T N (t): t > 0} then being defined by <T N (t) = exp (Ct N t), t > 0). With T N = <T N (x) and B N = (I - 
fl'N«> P N r+ + J 0 tr N ( s )P N Ar 4 ‘ds, the necessary conditions for convergence are argued using the 
well known Trotter-Kato semigroup approximation result (see (Kato, 1966)). The approximations 


* A A a 

to Q, Q and C,Qn, Qn and C N respectively, are typically taken to be Q N = P N Q, Q N = P N Q 
and 

C N = CP N ‘ 


Let { (pN} n N denote a basis for V N and set = (<p^ ,(p N ,..., <p£ ) T e x V N . 

j j =1 2 N j=i 

Adopting the convention that [L] denotes the matrix representation with respect to the basis 


n. 


{ tP N ) . N j for a linear operator L with domain and/or range in v N > we find that 
[F n ] = (R + [B n ] t 0 N [B n ]) 1 [B n ] t 0 n [T n ] and [F N ] = [T N ] ^[C^iR + [C N ] 


© n [C n ] t ) 1 where © N and 0 N are the unique, symmetric, nonnegative solutions to the 
n N xn N matrix algebraic Riccati equations 
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(3.10) 


0 N = [T n ] t (0N - 0 n [B n ](R + [B n ] t © n [B n ]) 1 [B n ] t e^)[T n ] 

+ M n [Q n ] 
and 

(3.11) 0 N = [T n ]( © N - © n [C n ] t (R + [C n ]0 n [C n ] t )' 1 [C n ]0 n )[T n ] t 

+ [Qn](M n )' ! . 

The matrix M N is the n N x n N Gramian matrix «X> N , (<X> N ) T > V . 

If we write w£j , k = (O 1 ^) 7 k with W* R 1 ^ , we obtain u* -[F N ] W Nk , 
k = 0,1,2,... with 


w H,..r ™ w h.> +[b ni u ;,, + ' f nh v - [c n ] k =°' 1 ' 2 '- 


W* =(M N )-’ <0", w„> v . 

N,0 u v 


The approximating optimal functional feedback control gains, f** = (f£ f£,..., f^) T e x V N 

j=i 

are given by f^ = [F n ](M n )* 1 O n and the approximating optimal functional observer gains 

fN=(fN %..., fN) T e x P V N by f = [F n ] T 0 N . Ifn N P N II weakly (strongly) 

j=l 

A A 

then f^ — » fj, i = 1,2,..., m weakly (strongly) in V. If II N P N —> 17 weakly (strongly) then 
f^ — > fj ,i = 1,2,..., p weakly (strongly) in V. If the injection VcH is compact, then in either 
case we have f|* -» fj . i = 1,2,..., m and fj . i = l,2,...,p strongly in H. 
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4. Examples and Numerical Result 


We consider the one dimensional heat equation 

n d^ w 

(4.1) (t,x) = a — - (t,x), 0 < x < l,t > 0, 

5t dx 

where a e R 1 , with a homogeneous Dirichlet boundary condition at x = 0, 

(4.2) w(t,0) =0, t > 0, 
and either Neumann, 

(4.3) (t,l) = v(t), t > 0, 

dx 

or Dirichlet, 

(4.4) w(t,l)=v(t), t>0, 

boundary control or input at x = 1 where v e 1^(0, <»). For output we take a temperature 
measurement 

(4.5) y(t) = w(t,C), t>0, 

at some fixed point t, e (0, 1). Initial conditions for these systems are assumed to be of the form 

(4.6) w(0, x) = w 0 (x), 0 < x < 1 

where w 0 e L 2 (0, 1). 

Although the two control systems above appear to be similar, they are, in fact, quite different and 
must be treated separately. We begin with the more straight forward of the two - the case of 
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Neumann boundary control. Let H = I^ (0, 1), V = H[ (0, 1) = {cp e FFfO, 1) : <p (0) = 0} and 
W = H 2 (0, 1)0 Hl( 0, 1). With H endowed with the usual L 2 inner product, V with the inner 
product < 9 , \j/ > v = Jq D 9 D\j/ and W with the inner product < 9 , y >w -/J D 2 9 D 2 \j/ we have the 
continuous and dense embeddings WcVcHcV'c W'. Define A e £(W, H), T e £(W, R 1 ) by 
A 9 = a 2 D 2 9 , r 9 = D 9 (l) and C 9 = 9 (C) respectively. With these definitions it is immediately clear 
that the boundary control system (4.1) - (4.3), (4.5), (4.6) is of the general form (2.1) - (2.4). The 
operator Cl: Dom (Cl) c H— >H is given by Cl 9 = a 2 D 2 9 for 9 e {9 e H^O, 1 ): 9 ( 0 ) = D 9 O)) = 

0} . It is densely defined, negative definite, self-adjoint and it is the infintesimal generator of a 
uniformly exponentially stable analytic semigroup { <T(t): t > 0 } of bounded, self-adjoint linear 
operators on H. In addition, { <T(t): t > 0} is also a uniformly exponentially stable, analytic 
semigroup of bounded, self-adjoint operators on V with generator £l given by £19 = Cl 9 for 9 e {9 
e H 3 (0, 1 ): 9(0) = D 9 (l) = D 2 9(0) = 0}. Choosing T+ e ^(R 1 , W) as ( T + u)(x) = xu for x e [0, 

1], we have fl\,( I 4- ) cV, R( T 4 ) c 71(A) and that conditions 1) -5) given in Section 2 are 
satisfied. For the optimal control and estimator problems, we take Q = ql, Q = q I, R = r and R = r 

A A 

where I is the identity operator on V, q, q >0 and r, r >0. The uniform exponential stability of the 
semigroup { S' (t): t > 0} on V implies that the algebraic Riccati equations (3.1) and (3.5) admit 

A 

unique, nonnegative self-adjoint solutions II and II respectively. The optimal control ( 3 . 2 ) takes the 
form 

1 

(4.7) u * = - J Df Dw* , k= 0 , 1 , 2 ,... 

0 

where the optimal functional feedback control gain f along with the optimal functional observer gain 
f are elements in H[ ( 0 , 1 ). 

We construct an approximation scheme using a linear spline based Ritz-Galerkin approach. For 
each N = 1,2,... let ( 9 ^}^_ Q denote the usual linear spline or "hat" functions defined on the 
interval [0,1] with respect to the uniform mesh {0, 1/N, 2/N,...,l } . We discard the element centered 
at x = 0, 9 ^, set V N = { 9 ^} l and choose P N to be the orthogonal projection of V onto V N 
with respect to the V inner product. It is clear that V N is an N dimensional subspace of V. 
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For (p e Dom (Cl), ICicpl H > a 2 l(pl v > a 2 l(pl H and therefore 0 e p(CO and Cl'*: H -» Dom(Ci) 

satisfies ICl'*<pl v < a' 2 l<pl H for <p e H. We define Ci N : V N -» V N as the inverse of the operator 

Cl' 1 = P N Cr J restricted to V N - That the operator Cl N is well defined follows from the fact that 
N 

”1 2 

< <Pn> v - -a' 2 l(p N l H for cp N e V N and that it is self-adjoint from <Ci N <p N ■ Vn > v = 

— a n^n^H 1 For (p N 6 Vjyj and = 9 estimate 

< ( ^n < Pn’ < Pn > v ~ < Vn^ n Vn>v = ' a ' 21 Vn 1 ^ 

< ^ICtVn'v^ ^IPnCtWn'J = -a 2 t Vn'v 

= ~ a2 I 9 nI v 


implies that Cl N is the infinitesimal generator of a Tjq semigroup { 0" N (t) : t > 0} of bounded, 

2 

self-adjoint linear operators on v N satisfying l<T N (t)l < e' a 1 , t>0. 

Let d N denote the interpolation operator from V onto V N- Then for (p e W, elementary 
approximation properties of linear in terpolatory spline functions (see (Schultz, 1971)) imply 

l(P N - I)(pl v < l(3 N - I)(pl v S — ID 2 <p| H 

Ntc 

and therefore, since W is dense in V, that P N — » I strongly on V as N— » oo . Also, it follows that Cl ^ 
= Pj^Cr 1 — » Cl' 1 strongly on V as N-» «» . If we define T N = <T n (t), then the Trotter-Kato 

ift 

approximation theorem (see (Pazy, 1983)) yields T N P N — » T (and therefore of course that T N P N — > 
T*) strongly on V as N— > °° where T = T* = 3" (x). 

Since <fi,(r + ) c V N , we define the approximating input operators B N by B N = (I - cT N (x))r + and 

A 

set Q n = ql, Q n = q I and C N = C. The strong convergence of P N to the identity and T N P N to T 
together with the finite dimensionality of the domain of B and the range of C are sufficient to 

A A 

conclude that Q N P N — > Q, Q N P N — » Q strongly on V and that B N — > B and C N P N — > C in norm as 
N — » 

The uniform exponential stability of the semigroups { <T N (t) : t > 0} implies 
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(4.8) IT* i v = l(T* )\ < r*. k = 0,1,2,... 


2 * 
with r = e ' a x < 1. Consequently the pairs (T N ,B N ) and (T N , C ) are uniformly exponentially 

A 

stabilizable and the pairs (T N ,Q N ) and (T N , Qn> are detectable, it follows therefore that there exist 
unique self-adjoint, nonnegative solutions II N and I1 N to the finite dimensional algebraic Riccati 
equations (3.7) and (3.8) respectively. The uniform exponential bound (4.8) with r < 1 on the 
approximating open-loop state transition operators t n imples that the zero control yields a uniform 

XV XV 

upper bound for II N and I1 N and the uniform exponential stability of S N = T N - B N F N and S N = T N 

A A A 

- F n C n . We can conclude therefore that n N Pn and n N P N converge strongly in V to I1 N and IT N 

XV A 

respectively and that P N P N and F n converge to F and F in norm as N — » «x>. The approximating 

XV A 

optimal functional feedback control and observer gains, f N and f N , converge respectively to f and f 
in the H 1 norm as N — » °°. 

In implementing the scheme outlined above, eigenvector decomposition of the associated 
Hamiltonian matrix was used to solve the matrix algebraic Riccati equations (3.10) and (3.1 1) (see 
Pappas, et. al., 1980). Matrix exponentials, where required, were also computed using 
eigenvalue/eigenvector decomposition. All calculations were carried out via codes written in Fortran 
and run on an IBM PC AT. We set a 2 = .1, q = q = r = ?= 1.0, £, = v2 /2 and x = .01 and obtained 

XV XV 

the functional gains plotted in Figs. 4.1 and 4.2 below. We plot f N and f N as well as Df N and Df N 
so as to exhibit the H 1 convergence and since it is Df (or Df N ) which actually appears as the feedback 
kernel in the optimal control law (4.7). 

We also simulated the operation of the closed-loop system with an approximating compensator. 
Using a 20 mode model for the infinite dimensional system and N = 12, we computed the closed- 
loop spectrum of the approximating compensator (i.e. the eigenvalues of the operator /6 N given by 
(3.9) with N = 12). These eigenvalues along with the first 20 open-loop eigenvalues (i.e. the first 
20 eigenvalues of the operator T = <T (x)) and the approximating closed-loop control and observer 
eigenvalues are tabulated in Table 4.1 below. Table 4.1 reveals that the last seven open-loop 
eigenvalues remain essentially unchanged in the closed-loop system-i.e. these modes are neither 
controlled nor observed by the finite dimensional compensator. Also, as one would expect, o(/6 N ) 
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consists essentially of the union of o(S N ) , a(S N ) and the eigenvalues corresponding to the 
uncontrolled/unobserved modes of the open-loop system. 

It is worth noting that the scheme we have outlined above for the Neumann boundary control 
problem is the same scheme that one would ordinarily use if the problem were formulated in the space 
H - i.e. if the output operator C was bounded on 1^(0, 1) (see Gibson and Rosen, 1986). This is 
possible primarily because the space V = H^(0,1) is the natural energy space for the underlying 
homogeneous or open-loop system. Consequently, the inherent self-adjointness and coercivity in the 
problem is preserved when it is formulated in the stronger space. In the case of Dirichlet boundary 
control, the situation is quite different. 

For the Dirichlet boundary control system (4.1), (4.2), (4.4) - (4.6), we choose the spaces H, V 
and W and their corresponding inner products to be the same as they were in the Neumann case. 

The operators A e £(W,H) and C e ^(V.R 1 ) also remain unchanged, however now we have 
r e ^(W.R 1 ) given by Tcp = cp(l). It then follows that the operator Cl : Dom(Cl) c H — » H is given 
by Cup = a 2 D 2 cp for cp e H 2 (0,1) fl H* (0,1). It is well known that Cl is densely defined, negative 
definite and self-adjoint and that it is the infinitesimal generator of the uniformly exponentially stable 
analytic semigroup { <T(t) : t > 0} of bounded, self-adjoint linear operators on H. However this 
time the operators <T (t) for t > 0 are neither self-adjoint nor a semigroup on V. Indeed, since 
ifi»(<r (t)) c H* (0,1) for all t > 0 and since H^(0,1) is a closed proper subspace of (0,1), <T (t) is 
not strongly continuous in the V-norm at t = 0. (The fact that our general framework requires 
IT+ = 1 and ^(P 1- ) c V precludes our choosing V to be H^(0,1).) On the other hand, 

{ <T(t) : t > 0} an analytic semigroup implies (see Pazy, 1983) that there exists a constant p > 0 for 
which ICl <T (t)l H < |lt~ 1 for t > 0 . Consequently, if we define T = <T (t), then it follows that 
T e £(V) and moreover, that 


IT k (pl 2 = - <Cl O' (kx)<p,<T (kx)<p> < I CiiT(kt)(pl H I tT (kx)<p| H 


pe 


2 

-a kt 


kp? H 


pe 


2 

- a kx 




kx kx 

for k = 1,2,... and <p e V. We have therefore 
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(4.9) IT k l v = l(T*) k l v < MrK k = 0,1,2,... 


where M > 0 and r < 1. 

We again choose e ^(R^W) as ( r f u)(x) = xu for x e [0,1]. Then lfk(r + ) c rj(A) and we 
have reformulated the boundary control system (4.1), (4.2), (4.4) - (4.6) in the general form of (2.1) 
- (2.4) and conditions 1) - 5) are satisfied. 

We formulate the optimal control problem with the performance index 


-2 


j(u)= Zj q <w k ’ w k > H +ru k 

k=0 

where q > 0 and r > 0. That is, we take Q to be the bounded, self-adjoint nonnegative operator on 
Hl( 0,1) given by (Qcp)(x) = q Jq J y <p(z)dzdy and R to be r. For the estimator problem we set 

A A /V A 

Q = q I and R = r with q > 0 and r > 0. 

The uniform exponential bound (4.9 ) implies the existence of unique, nonnegative, self-adjoint 

A 

solutions n and IT to the algebraic Riccati equations (3.1) and (3.5). The optimal control is again of 

~ 1 

the form (4.7) with the optimal functional gains f and f in H L 

The fact that { iT(t) : t > 0} is not a semigroup on V precludes the use of a semigroup - theoretic 

approach to approximation. We therefore employ modal subspaces and approximate the open-loop 

state transition operator T directly as a bounded linear operator on V. 

For each N = 1,2,... let V N = span { (pj } q where for x e [0,1], <p 0 (x) = x and <pj(x) = sinjrcx, 

j = 1,2,...,N. Let p N denote the orthogonal projection of H = L 2 (0,1) onto span {(p;} N 

J j =i 

and let P N denote the orthogonal projection of V onto v N - Using the fact that V = H*(0,1) ® <Pq, it 
is not difficult to see that P N tp = <p(l)<p 0 + Pn( 9 “ 9(1)%) for <p e V and hence, via elementary 
properties of Fourier series (see Tolstov, 1962), that l(P N - I)(p|y = l(p N - I)(cp - <p(l)%)| v — > 0 

as N ©o for each <p e V. 

N 

If for V N = X 9j e V N we define T N e ^(V N ) by 
i=o 

N j 2 2 

t n V n = P n t ¥ n = P n ^(x)V n = ]£ { ^- VN + V N )e' “ J n2 \ r 

j=i 
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then = P N T*, iT k I = l(T^ k l y < Mr k , k = 0,1,2,... with M > 0 and r < 1 independent of N, 


l(T N P N ■ T)cply < KPjsjTPjsf ■ PjsjTl^pIv + KP n ■ I)Tcply 
< Mrl(P N - I)cpl v + l(P N - I)Tcp| v -> 0 


as N oo for cp e V. Similarly, r p N -4 t* strongly on V as N °o. 

The approximating input, output and state penalization operators B N , C N , Q N and Q N take the 
form 


<5/1 d 222 

B N U =(I-T N )r + u =cp Q u+ — — e' a 1 n T (p.u , 


j=i JJt 


C N ~ C,Q n = qP N Q and q n = qi- Reasoning as we did in the Neumann case, the approximating 
algebraic Riccati equations (3.7) and (3.8) admit unique, nonnegative, self-adjoint solutions n N and 

A A A A A 

n N respectively, n N P N —> II and n N P N - 4 II strongly on V and F N P N — > F and F N -4 F in norm 

A 

as N -4 00 . The approximating functional feedback control and observer gains f N and f N converge 

A 

to f and f respectively, strongly in H 1 as N -4 <*> . 

With a 2 = 1.0, q = q = r = 1.0, r = 5.0, £ = V2/2 and x = .01 and the scheme outlined above we 
obtained the approximating optimal functional feedback control and observer gains plotted in Figs. 

4.3 and 4.4 below. The first 12 open-loop and the approximating closed-loop control and observer 
eigenvalues for N = 12 are tabulated in Table 4.2. 

Table 4.2 reveals an interlacing of the closed-loop control and open-loop eigenvalues. That is, 
the closed-loop control eigenvalues (i.e. the elements in the spectrum of S) are alternately more and 
less stable than the corresponding open-loop eigenvalues. We also have observed this phenomenon 
in other numerical studies we are carrying out involving LQG boundary control for flexible 
structures. In additon, in the Dirichlet boundary control system discussed above, if Q is chosen as 
the identity operator on V = H^(0,1), virtually all of the closed-loop control eigenvalues are less stable 
than the corresponding open-loop eigenvalues. It is clear that this non-standard behavior results from 
the presence of the one dimensional subspace represented by U\»(r + ). Indeed, the behavior of the 
closed-loop spectrum in the case of Neumann boundary control is as would be expected. We feel that 
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what we are seeing can most likely be explained via infinite dimensional analogs of existing results 
relating the asymptotic properties of the closed-loop spectrum of a linear regulator and the zeros of 
the corresponding open-loop transfer function (see Kwakemaak and Sivan, 1972 and Harvey and 
Stein, 1978). However, as of yet, we have been unable to establish this conjecture satisfactorily 
and we consider it to be beyond the scope of this paper, which is primarily concerned with 
approximation. We leave it as an interesting open question worthy of further investigation. 


5. Concludine Remarks 

We have developed a framework for the finite dimensional approximation of optimal discrete-time 
LQG compensators for distributed parameter systems with boundary input and unbounded 
measurement. Our theory applies to the class of boundary control problems which can be formulated 
in a state space in which both the discrete-time input and output operators are continuous. We have 
used a functional analytic treatment to develop a convergence theory and have demonstrated the 
feasibility of our approach via examples involving either the Neumann or Dirichlet boundary control of 
a one dimensional heat equation with point measurement of temperature. We have shown that while 
both problems outwardly appear to be quite similar, they in fact require very different approaches to 
approximation. Also in the Dirichlet case the observed behavior of the resulting closed-loop spectrum 
is, in some ways unexpected and its explanation remains open. 

Finally, we have been looking at the application of our schemes to LQG problems for flexible 
structures with boundary inputs and unbounded measurement and systems with control and/or 
observations delays. We have been considering vibration suppression for cantilevered beams via 
shear or moment inputs at the free end and pointwise observation of strain or acceleration. These 
studies are currently underway with the results to be reported elsewhere. 

Acknowledgment: The authors would like to gratefully acknowledge Mr. Milton Lie of the 
Department of Mathematics at the University of Southern California for his assistance in carrying out 
the computations reported on in this paper. 
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Open-Loop 


.9975 

.9780 

.9402 

.8861 

.8188 

.7419 

.6590 

.5740 


.4104 

.3368 

.2711 

.2139 

.1655 

.1255 

.0934 

.0681 

.0482 

.0341 

.0235 


.9968 .9968 

.9780 

.9769 .9768 

.9408 

.9371 .9371 

.8872 

.8778 .8775 

.8194 

.7982 .7985 

.7414 

.7026 .7019 

.6573 

.5960 .5921 

.5718 

.4891 .4769 

.4875 

.4433 

.4041 

.3675 .3675 

.3341 

.2772 .2763 

.2705 

.2145 .2129 

.2134 

.1811 .1811 

.1663 
.1260 
.0933 
.0677 
.0483 
.0340 
.0236 


Neumann boundary control; simulation results 


Table 4.1 
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Open-Loop 

a(Si 2 ) 

a(Si 2 ) 

.90601806 

.90569591 

.78573771 

.6 7382545 

.68243047 

.57981918 

.41136911 

.40961171 

.40082268 

.20615299 

.20758391 

.20936323 

.08480497 

.08447005 

.08636884 

.02863695 

.02873534 

.02892353 

.00793790 

.00791793 

.00792193 

.00180617 

.00180978 

.00178763 

.00033753 

.00033682 

.00033414 

.00005172 

.00005179 

.00005162 

.00000651 

.00000650 

.00000654 

.00000067 

.00000067 

.00000068 


.00000000 

.00000000 


Dirichlet boundary control; open and closed-loop spectrum 

Table 4.2 
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